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Abstract 

A number of fields, including the study of genome rearrangements 
and the design of interconnection networks, deal with the connected 
problems of sorting permutations in "as few moves as possible" , using 
a given set of allowed operations, or computing the number of moves 
the sorting process requires, often referred to as the distance of the 
permutation. These operations often act on just one or two segments 
of the permutation, e.g. by reversing one segment or exchanging two 
segments. The cycle graph of the permutation to sort is a fundamental 
tool in the theory of genome rearrangements, and has proved useful 
in settling the complexity of many variants of the above problems. 
In this paper, we present an algebraic reinterpretation of the cycle 
graph of a permutation ir as an even permutation W, and show how to 
reformulate our sorting problems in terms of particular factorisations 
of the latter permutation. Using our framework, we recover known 
results in a simple and unified way, and obtain a new lower bound on 
the prefix transposition distance (where a prefix transposition displaces 
the initial segment of a permutation), which is shown to outperform 
previous results. Moreover, we use our approach to improve the best 
known lower bound on the prefix transposition diameter from 2n/3 to 
|_3n / 4 J , and investigate a few relations between some statistics on tt 
and W. 



*A significant portion of this work previously appeared in the Proceedings of the Six- 
teenth Annual European Symposium on Algorithms (ESA) (22|. 
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1 Introduction 



Given a set S of allowed operations and two permutations 7r and a of {1, 2, . . . , n}, 
we study the related problems of computing, on the one hand, a sequence of 
elements of 5* of minimum length that transforms n into a, and on the other 
hand, computing the length of such a sequence, referred to as the distance 
between ti and a. The operations in 5* usually yield an edit distance ds{-,-) 
with the property that ds(n,cr) = dsi^cr^ 1 o ir, l) for any two permutations 
7r and a of the same set, where i is the identity 'permutation (12 • • ■ n). 
This property allows us to restrict our attention to sorting permutations us- 
ing a minimum number of operations from 5*, or to computing the distance 
of a given permutation to the identity permutation rather than to another 
arbitrary permutation. Two areas in which these questions have applications 
are the fields of genome rearrangements and interconnection network design, 
which we briefly review below. 

In genome rearrangements (see Fertin et al. 12[ for a survey), the per- 
mutation to sort represents an ordering of genes in a given genome, and 
the allowed operations model mutations that are known to actually occur 



in evolution. Rearrangements studied in that context include reversals [19 



which reverse a segment of the permutation, block-interchanges [8(, which ex- 
change two not necessarily contiguous segments, and transpositions [i[ , which 
displace a block of contiguous elements. Those seemingly easy problems 
turn out to be more challenging than they might appear at first: although 
a polynomial-time algorithm is known for sorting by block-interchanges or 
computing the associated distance jsj, the same problems were^shown to be 
NP-hard for reversals 0], and more recently for transpositions 



In interconnection network design (see Lakshmivarahan et al. |23| for a 
thorough survey), permutations stand e.g. for processors, or other devices 
to be connected, and form the vertex set of a graph whose edges corre- 
spond to physical connections between two devices. One wants to build a 
grap h with small degree and sm all diameter, among other desirable proper- 
ties. lAkers and Krishnamurthyr s landmark paper |l| proposed the idea of 
choosing a set S that generates all permutations of {1, 2, ... , n}, and to use 
the corresponding Cayley graph, whose vertex set is the set of all permuta- 
tions and whose edges connect any two permutations that can be obtained 
from one another by applying a transformation from 5*, as an interconnection 
network. In that setting, sorting algorithms for permutations correspond to 
routing algorithms for the corresponding networks, since a sequence of ele- 
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ments of S transforming it into a corresponds to a path of the same length 
in the network. Two kinds of operations that received a lot of attention in 
that context are prefix reversals [l5|], which reverse the initial segment of the 
permutation, and prefix exchanges [2j|, which swap the first element of the 
permutation with another element. Those operations gave birth to the pan- 
cake network and star graph topologies, respectively, which are extensively 
studied models in that field. We also mention prefix transpositions, which 
displace the initial segment of the permutation, and were introduced by Dias 
and Meidanis jsj in the context of genome rearrangements in the hope that 
their study would shed light and give insight on the challenging problem of 
sorting by transpositions. Those more restricted versions of operations stud- 
ied in the context of genome rearrangements do not lead to problems simpler 
than their unrestricted counterparts: the sorting and distance computation 
problems related to prefix exchanges can be solved in polynomial time 
but the complexity of those problems in the case of prefix transpositions is 
open, and the problem of sorting by prefix reversals has only recently been 
showed to be NP-hard jij], more than thirty years after the first works on the 



subject |15|, [16 



The cycle graph of a permutation is a ubiquitous structure in the field of 
genome rearrangements, and has proved useful in resolving many questions 
related to the problems discussed in the above paragraphs. In this paper, we 
present a new way of encoding the cycle graph of a permutation tt as an even 



permutation W, inspired by a previous work of ours [10|, and show how to 
reformulate any sorting problem of the form described above in terms of par- 
ticular factorisations of the latter permutation. We first illustrate the power 
of our framework by recovering known lower bounds on the block-interchange 
and transposition distances in a simple and unified way, and then use it to 
prove a new lower bound on the prefix transposition distance. We prove that 
our lower bound always outperforms that obtained by Dias and Meidanis [of , 
and show experimentally that it is a significant improvement over both that 
result and the only other known lower bound proved by Chitturi and Sudbor- 
ough [7]. We then use this new result to improve the previously best known 
lower bound on the maximal value of the prefix transposition distance from 
2n/3 to [3n/4j. Finally, we examine some further properties of the model, 
and establish connections between statistics on ti and if. 
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2 Notation and definitions 



2.1 Permutations and conjugacy classes 

Let us start with a quick reminder of basic notions on permutations (for 
details, see e.g. Wielandt H)- 

Definition 2.1. A permutation of a set Q is a bijective application of Q onto 
itself. 

It is convenient to set Q = {1,2, ... ,n}, and we will follow this conven- 
tion here, although we will also sometimes use the set {0, 1, 2, ... , n}. The 
symmetric group S n is the set of all permutations of a set of n elements, 
together with the usual function composition o, applied from right to left. 
Permutations are denoted by lower case Greek letters, and we will follow the 
convention of shortening the traditional two-row notation 

/ 1 2 ••• n 

7T = ( 

\ 7Ti 7T 2 ■ • • 7T n 

by keeping only the second row, i.e. n = (n 1 n 2 ■ ■ ■ ir n ), where 7Tj = 7r(z). 

Definition 2.2. The graph T(tt) of the permutation 7r in S n is the directed 
graph with ordered vertex set (tti, 7r 2 , . . . , n n ) and arc set {(i, j) | 7Tj = j, 1 < 
i < n}. 

Definition 12.11 implies that r(7r) decomposes in a single way into disjoint 
cycles (up to the ordering of cycles and of elements within each cycle), leading 
to another notation for ir based on its disjoint cycle decomposition. For 
instance, when n = (4 16257 3), the disjoint cycle notation is n = 
(1, 4, 2)(3, 6, 7)(5) (notice the parentheses and the commas). 

Definition 2.3. The length of a cycle in a graph is the number of vertices 
it contains, and a k-cycle is a cycle of length k. 

The number of cycles in a graph G will be denoted by c{G), and the 
number of cycles of length k will be denoted by CkiG). We will also distinguish 
between cycles of odd (resp. even) length, denoting the number of such cycles 
in G using c dd(G) (resp. c even (G)). It is common practice to omit 1-cycles 
in the cycle decomposition of (the graph of) a permutation, and to call that 
permutation a fc-cycle if the resulting decomposition consists of a single cycle 
of length k > 1. Cycles of length 1 in the disjoint cycle decomposition of a 
permutation are referred to as fixed points. 
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Definition 2.4. A permutation n is even if the number of even cycles in 
r(7r) is even or, equivalently, if it can be expressed as a product of an even 
number of 2-cycles. 

The alternating group A n is the subgroup of S n formed by the set of all 
even permutations, together with o. The following notion will be central to 
this work. 

Definition 2.5. The conjugate of a permutation 7r by a permutation a, 
both in S n , is the permutation ix u = a o n o c -1 , and can be obtained by 
replacing every element i in the disjoint cycle decomposition of n with o~{. 
All permutations in S n that have the same disjoint cycle decomposition form 
a conjugacy class (of S n ). 

2.2 Generating sets and edit distances 

We are interested in distances between permutations based on operations 
that can themselves be modelled as permutations. More formally, given a 
subset S of S n and two permutations 7r and a in S n , we have two goals: 

1. to find a sequence of elements si, S2, ■ ■ ■ , St from S whose cardinality 
is minimum and whose product transforms n into a (or conversely, a 
into 7r) : 

7T O Si O S 2 O • • • O S t = O. 

2. to find the length of such a sequence, called the S distance between 7r 
and a. Distances whose definition is based on a set of allowed opera- 
tions as described above are often referred to as edit distances. 

Note that S must be symmetric, i.e. s e S if and only if s _1 e 5, for 
the corresponding distance to satisfy the symmetry axiom. An immediate 
corollary of this property is that for any n in S n , we have d(ir, l) = din^ 1 , t). 
For any two permutations of the same set to be a finite distance apart, S 
must also satisfy the following property. 

Definition 2.6. A set S C S n is said to generate S n , or to be a generating 
set of S'n, if every element of S n can be expressed as the product of a finite 
number of elements of S. We call the elements of S generators of S n . 

Moreover, all generating sets we will consider in this paper yield distances 
that satisfy the following property. 
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Definition 2.7. A distance d on S n is left-invariant if for all tt, a, r in S n , 

we have: d(ir, a) = d(r o tt, r o a). 

Intuitively, left-invariance models the fact that, given any two permuta- 
tions 7r and a to be transformed into one another, we can rename the elements 
of either permutation as we wish without changing the value of the distance 
between both permutations, as long as we renumber the elements of the other 
permutation accordingly. Since we will most of the time be considering the 
distance between a permutation tt and the identity permutation l, we will 
often abbreviate d(ir, l) to d{rr). 

It can be easily seen that both problems mentioned at the beginning 
of this section can be reformulated in terms of finding a minimum-length 
factorisation of tt that consists only of elements of S, since 

7TOSlOS2 0- -- OS ( = t<^7r = s^ 1 o o • • • o s^ 1 

and S is symmetric. Finally, another parameter of interest in the study of 
those distances is the largest value they can reach. 

Definition 2.8. The diameter of a set U under a distance d is max Sjte j/ d(s, t). 



2.3 Genome rearrangements and the cycle graph 

We recall here a few operations that are commonly used in the fields of 
genome rearrangements and interconnection network design to build gener- 
ating sets of S n . 

Definition 2.9. [8j The block-interchange (3(i,j,h,l) with l<i<j<k< 
I < n + 1 is the permutation that exchanges the closed intervals determined 
respectively by i and j — 1 and by k and I — 1: 



1 ■ 


• i- 1 


% ■ ■ 


■ 3-1 




■ jfc-1 


k ■ 


■■ I - 1 


ll + l ■ 


■ n 


1 ■ 


• i- 1 


k ■ 


■■ I- 1 


3 3 + 1 ■ 


• k- 1 


i ■ 


■3-1 


ll + l ■ 


■ n 



Two particular cases of block-interchanges are of interest: 



1. when j = k, the resulting operation exchanges two adjacent intervals, 
and is called a transposition [3|, denoted by r(i,j,l); 
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2. when j = i + 1 and I = k + 1, the resulting operation swaps two not 
necessarily adjacent elements in respective positions i and k, and is 
called an exchange, denoted by e(i,k). 

We use the notation bid(7r), td(ir) and exc(ir) for the block-interchange 
distance, the transposition distance, and the exchange distance of 7r, respec- 
tively. The operations we described above can be further restricted by setting 
i — 1 in their definition, thereby transforming them into so-called "prefix 
rearrangements". The corresponding "prefix distances" are defined in an 
analogous manner, with the additional restriction that all operations must 
act on the initial segment of the permutation. We denote ptdiff) and pexc(ir) 
the prefix transposition distance and prefix exchange distance of ir, respec- 
tively. While sorting by transpositions is NP-hard |4( and the computational 
complexity of sorting by prefix transpositions is unknown, polynomial-time 



algorithms exist for sorting by block- interchanges |8j , exchanges [18[ or prefix 
exchanges jij, as well as formulas for computing the associated distances. 

We will have more to say about sorting by transpositions and sorting by 
block-interchanges in Section HI where we will give simple proofs of lower 
bounds on the two corresponding distances, as well as about sorting by 
prefix transpositions in Section El where we will prove new and improved 
lower bounds on the associated distance and diameter. Meanwhile, we con- 
clude this section with the following traditional tool introduced by Bafna and 
Pevzner j^] , which has proved most useful in the study of genome rearrange- 
ments. 

Definition 2.10. The cycle graph of a permutation tt in S n is the bicoloured 
directed graph G(n), whose vertex set (7r = 0, 7Ti, . . . , 7r n ) is ordered by 
positions, and whose arc set consists of: 

• black arcs {(7Tj,7Tj_i) | 1 < i < n} U {(7To,7r n )}; 

• grey arcs {{^i, 7Tj + 1) | < i < n} U {(n, 0)}. 

The arc set of G(n) decomposes in a single way into arc-disjoint alter- 
nating cycles, i.e. cycles that alternate black and grey arcs. The length of 
an alternating cycle in G(ir) is the number of black arcs it contains, and a 
k-cycle in G{ir) is an alternating cycle of length k (note that this differs from 
Definition 12. 3j) . Figure [1] shows an example of a cycle graph, together with 
its decomposition into a 5-cycle and a 3-cycle. 
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Figure 1: (a) The cycle graph of (4 1 6 2 5 7 3); (6), (c) its decomposition 
into two alternating cycles. 

3 A general lower bounding technique 

We now present a framework for obtaining lower bounds on edit distances 
between permutations in a simple and unified way. To that end, we adapt a 
bijection previously introduced by Doignon and Labarre (lo| : 




/ : 5 n A n+l : 7r H> 7r = (0, 1, 2, . . . , n) o (0, vr n , vr n _i, . . . , vri), (1) 

which in particular maps t onto I = (0 12 • • • n). That mapping allows 
us to encode the structure of a cycle graph G(ir) using an even permutation 
7f in an intuitive way, which corresponds to decomposing the cycle graph 
into the product of two "monochromatic cycles" , namely, the cycle made of 
all black arcs (i.e. (0, n n , 7T n _i, . . . , 7Ti)) and the cycle made of all grey arcs 
(i.e. (0, 1, 2, . . . , n)). The construction is perhaps best understood using an 
example: let tt = (4 1 6 2 5 73), whose cycle graph is depicted in Figure [Qa). 



and the two disjoint cycles of Tf correspond to the two alternating cycles 
of G(tt), whose elements they list in the order they are encountered (up to 
rotation); indeed: 

1. the first cycle of G(ir) (Figure []](&)) starts with 0, then visits 4 after 
following a black-grey path (i.e. a black arc followed by a grey arc), 




Then 



7T = 



(0, 1, 2, 3, 4, 5, 6, 7) o (0, 3, 7, 5, 2, 6, 1, 4) = (0, 4, 1, 5, 



3)(2,7,6), 
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then visits 1 after following a black-grey path, and in the same way 
visits 5 and 3 before coming back to 0, which corresponds to the first 
cycle of 7f ; 

2. the second cycle of G(tt) (Figure Q^c)) starts with 2, then visits 7 after 
following a black-grey path, and in the same way visits 6 before coming 
back to 2, which corresponds to the second cycle of W. 

Note that the order in which we decide to follow arcs (first a black 
arc and then a grey arc) is given by the order in which the two cycles 
are multiplied. An alternative definition^ of W could therefore have been 
(0, 7r n , 7r n _i, . . . , 7Ti) o (0, 1, 2, . . . , n), which can be seen to be equivalent to 
our definition when conjugated by (0, n, n— 1, . . . , 1), and whose cycles are in- 
terpreted exactly as above, with the modification that grey arcs are followed 
first. Consequently, speaking about cycles of W, of r(7f) or of G(ir) is equiva- 
lent. We will now demonstrate how /(•) can be used to obtain results on the 
sorting and distance computation problems we discussed in Section 12.21 The 
following lemma expresses how the action of any rearrangement operation 
a on tt is translated on Jr. We will find it convenient to identify permu- 
tations in S n with their extended versions in S n+ i (i.e. we identify tt with 
(0 TV i 7T2 • ■ ■ 7T n )). This allows us to express any permutation tt in S n as 
follows: 

TT = (0,7T n ,7r n _i, . . . ,7Ti) OTTO (0, 1, 2, . . . , Tl) . (2) 

Lemma 3.1. For all tt, a in S n , we have tt o a = Jfo a 71 . 
Proof. By definition, we have: 

TToa = (0,l,2,...,n) O (0, (TTOa) n , (7T O Cr) n _i, . . . , (TT O (t)i) 

= (0, 1, 2, . . . , n) o 7r o (0, a n , er n _i, . . . , 0i) o tt' 1 

= (0, 1, 2, . . . , n) O (0, TT n , 7T„_i, . . ■ , 7Ti) O 7T o (0, 1, 2, . . . , n) 

o(0, a n , a n _i, . . . , oi) o tt' 1 (using Equation [2]) 

□ 

We are now ready to prove our main result. 



lr This is actually the definition we used in the conference version of this paper 22 1. 
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Theorem 3.1. Let S be a subset of S n whose elements are mapped by /(•) 
onto S' C A n+1 . Moreover, let be the union of the conjugacy classes (of 
S n+ \) that intersect with S ; ; then for any ir in S n , any factorisation of n into 
t elements of S yields a factorisation ofW into t elements ofio. 

Proof. Induction on t. The base case is tc G S, and clearly W G S' C . For 
the induction, let n = g t o g t _% o ■ ■ ■ o gi, where gi G S for 1 < i < t, and let 
a = g t -\ o ■ ■ • o g 2 o gi ; by Lemma I3.1[ we have: 



= 3t° 0t-i ° ■ ■ ■ ° g2 ° gi = gt ° & = g t o a* 1 . 

By induction, a = o o • • ■ o g[ , where g\ G ^ for 1 < i < t — 1; 
therefore: 

= ft° St-i ° g't-2 ° ■ ■ ■ ° 9i ° 9t X 

= gt ° g't-i ° £T^ ° ff* ° £t-2 ° ^ °fft --- ft lo 9t°g| 9t\ » 

and hi,h 2 , ■ ■ ■ , ht-i G ^, which completes the proof. □ 

We will use Theorem 13.11 in the next two sections to prove lower bounds 
on several edit distances between permutations. 



4 Recovering previous results 

We illustrate how to use Theorem l3.1l to recover two previously known results 
on bid and td. The general idea is as follows: as we explained in Section \272\ if 

5 is symmetric, then any sorting sequence of length t for 7r made of elements 
of S yields a factorisation of 7r into the product of t elements of S, which can 
in turn be converted, as in the proof of Theorem 13 .11 into a factorisation of W 
into the product of t elements of S' C % \ Therefore, the length of a shortest 
factorisation of 7f into the product of elements of ^ is a lower bound on the 
length of a factorisation of n into the product of elements of S, and we can 
obtain a lower bound on the distance of interest by: 

1. characterising the set of images of the elements in S by /(•), and 

2. computing the distance of W with respect to < rf. 
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Let us now show how we can obtain a lower bound on the block-interchange 
distance. We start by characterising the image of a block-interchange by our 
mapping. 

Lemma 4.1. For any block-interchange (3(i,j,k,l) in S n , we have 



P(i,j,kJ) = UJ) ° (i,k). 
Proof. Equation [1] and Definition 12.91 yield 

(0, 1,2, ...,n) o (0,n,n- 1, j - 1, j -2, ...,i,k- l,k- 2, 
j,l-l,l-2,...,k,i-l,i-2,...,l) 
= (0)(n)(n - 1) • • • (I + m,M ~ 1) • • • (i + k)(i - 1) • • • (1) 



□ 



Note that (j, I) and (i, k) might not be disjoint, since Definition 12.91 allows 
for j = k (hence the use of o in the expression of (3(i,j, k, /)). We can now 
recover a known lower bound on the block-interchange distance, which is 
actually the exact distance as shown by Christie 

«+i-c(r(7f)) 



Theorem 4.1. Ja/ For all it in S n , we have bid(7v) > 

Proof. By Theorem 13.11 and Lemma 14.11 a lower bound on bid(7r) is given 
by the length of a minimum factorisation of W into the product of pairs of 
exchanges. Since this length equals {n+ 1 — c(r(?f)))/2 (see e.g. Jerrum [181] ). 
the proof follows. □ 

Let us now characterise the image of a transposition by our mapping. 
Lemma 4.2. For any transposition r(i,j,l), we have 



r(i,j,l) = 

Proof. As noted in Section [2T2l we have r(i,j,l) = f3(i,j,j,l); Lemma H~T1 
yields 

□ 
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We can now recover the following known lower bound on the transposition 
distance. Recall that c ot y(r(7?)) denotes the number of odd cycles in r(?f). 

Theorem 4.2. [5] For all n in S n , we have td(ir) > n+1 ~ Co 2 dd(r(¥)) . 

Proof. By Theorem 13.11 and Lemma I4.2[ a lower bound on td{ir) is given 
by the length of a minimum factorisation of W into the product of 3-cycles. 
Since this length equals (n + 1 — c odd (T(li)))/2 (see e.g. Jerrum 0), the 
proof follows. □ 



New results on the prefix transposition dis- 
tance 



Dias and Meidanis [9] initiated the study of sorting by prefix transpositions, 
and derived a lower bound on the corresponding distance using the following 
concepts. 

Definition 5.1. Given a permutation tt in S n , build the permutation tt = 
(0 TV i • • ■ 7r n n + 1); a pair (7^, 7r i+1 ) with < i < n is a prefix transposition 
breakpoint if ^ ^ + 1 or if i = 0, and an adjacency otherwise. 

The number of prefix transposition breakpoints of tt is denoted by ptbfa). 
Noting that a prefix transposition can create at most two adjacencies and 
that i is the only permutation with one prefix transposition breakpoint, 
Dias and Meidanisl obtained the following lower bound. 



Lemma 5.1. |1/ For any it in S n : 

ptd(n) > 1 



(3) 



Chitturi and Sudborough 0| later obtained another lower bound on the 
prefix transposition distance. They used the following concepts, based on 
permutations of {0, 1, 2, . . . , n — 1} rather than {1, 2, . . . , n}. 

Definition 5.2. For a permutation 7T of {0, 1, 2, . . . , n — 1}, an ordered pair 
(7Tj, 7Tj + i) is an anti- adjacency if = 7^ — 1 (mod n). A strip in a permu- 
tation it is a maximal interval of 7r that contains only adjacencies, and a clan 
is a maximal interval of tt that contains only anti-adjacencies. 
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Chitturi and Sudboroughl proved the following lower bound. 



Lemma 5.2. For any permutation it o/{0, 1, 2, . . . , n— 1}, letTijx) denote 
the set of all clans of it of length at least 3 7 and s(n) denote the number of 
strips of n. Then 

Mt) > — 2 3 • (4) 

We will prove a new lower bound on the prefix transp osition distance 
(Theorem 15.21 page [TJJ, using our model and Akers et al. 's results on 
computing the prefix exchange distance: 

Theorem 5.1. 0/ For any ir in S n , we have 



pexc(n) = n + c{T{ti)) — 2ci(T(7r)) 



i/7Ti = l, 

2 otherwise, 



where Ci(T(n)) denotes the number of 1-cycles in T(n), or equivalently the 
number of fixed points of n. 



5.1 An improved lower bound 

Using our theory, we prove a new lower bound on ptd(ir) and show that 
it always outperforms ([3]). We will find it convenient to express ptb{ix) as 
follows. 

Lemma 5.3. For any tc in S n , we have 

p«w=«+i-c 1 (r W)+ {j f2;i 

Proof. The formula results from the observation that, among the n + 1 pairs 
of adjacent elements in n, each adjacency in 7r gives rise to a 1-cycle in r(7f), 
and from the fact that if ti\ = 1, then we counted the 1-cycle that corresponds 
to (0, 1) as an adjacency, which is contrary to Definition 15.11 and which we 
correct by adding 1. □ 

As explained in Section HJ we can obtain a lower bound on the prefix 
transposition distance by characterising the image of a prefix transposition 
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by /(•) and computing the associated distance. We already know that trans- 
positions are mapped onto 3-cycles (see Lemma [4.21 page ITTj) ; in the case of 
prefix transpositions, it is easily seen that these 3-cycles will always contain 
element 0. Therefore, we need to be able to compute the length of a mini- 
mum factorisation of tt in S n into a product of 3-cycles, where each 3-cycle 
in the factorisation is further required to contain the first element. Let us 
denote the corresponding distance dl(n); the following result shows how to 
compute it. 

Lemma 5.4. For any tt in A n , we have 

n + c(T(n)) f if <k x = 1, 



dm = o - <*( r fr)) - 



1 otherwise. 



Proof. Given a minimum factorisation of length i of an even permutation tt 
into prefix exchanges, we can construct a sequence of £/2 3-cycles by noting 
that (l,j) o = Therefore dl(ir) < £/2. On the other hand, 

assume there exists a shorter sequence of 3-cycles acting on the first element 
whose product is tt; then one can split each of these 3-cycles into two pre- 
fix exchanges using the relation above and find a shorter expression for tt 
as a product of prefix exchanges, a contradiction. The result follows from 
Theorem 15.11 □ 

As a corollary, we obtain the following new lower bound on the prefix 
transposition distance. 

Theorem 5.2. For any tt in S n , we have 

^)> " ti+ 3 ^" -,m)-(; iri (5 ) 

Proof. Follows immediately from Theorem 13.11 and Lemma 15.41 □ 

An immediate question is how tight this new lower bound actually is. 
We will answer this question experimentally in Section |6l where we will see 
that many more permutations are tight with respect to our new result than 
with respect to the previously known lower bounds. We will in the meantime 
conclude this section by proving that our lower bound always outperforms 
Dias and Meidanis' (given by Lemma [5. ip . 



Theorem 5.3. For all tt in S n , the value of lower bound (jSJ) always exceeds 
that of lower bound . 
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Proof. Assume tt ^ i (otherwise the result trivially holds); this implies that 
r(7f) has at least one cycle of length at least 2, which means that c(T(W)) — 
ci(r(7f)) > 1. There are two cases to prove: if 7Ti = 1, then lower bound (j3J) 
becomes 



(n + l- Cl (r(7f)) + l)-l 



n + l- Cl (r(7f))' 



and lower bound (jSJ) satisfies 

71 + 1 + C(r(7f)) -2Ci(r(7f)) 7l + 2- Ci(r(7f)) 



n + 1 - ci(r(7f)) 



On the other hand, if ~K\ ^ 1, then lower bound (|3]) becomes 



(n + l- Cl (r(7f)))-l 



n- Cl (r(7f)) 



and Definition l2.4l implies that for any n in A n , we have n = c(r(7r)) (mod 2). 
Lower bound ([5]) becomes 



n + 1 + c(r(vf)) 



Cl (r(7f))-i 



> 



n+l + c(r(7f))-2 Cl (r(7f))-2 



"n - ci(r(7r))' 



□ 



5.2 A tighter lower bound on the prefix transposition 
diameter 



Dias and Meidanis [9j observed that the prefix transposition diameter lies 
between n/2 and n — 1, and conjectured that it is equal ton - |_^J . Chitturi 
and Sudborough 0] then improved those bounds to 2n/3 and n — log 8 n, 
respectively. Using our new lower bound, we further improve the lower bound 
on the prefix transposition diameter. We prove our result in a constructive 
way, by building families of permutations whose prefix transposition distance 
is at least [3^/4] . Figure [21 which follows our result, shows examples of such 
permutations. The proof uses permutations from the following class, which 



has proved useful in the analysis of several other rearrangement problems [12 
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Definition 5.3. A permutation ir in S n is a 2-permutation if all cycles in ir 
have length 2. 

Note that the above definition requires n = 3 (mod 4): indeed, n + 1 
must be even in order to obtain a partition of the elements of if into pairs, 
and (n + l)/2 is also even by the definition of W. 

Theorem 5.4. For all n, the prefix transposition diameter of S n is at least 
[3n/4j . 

Proof. If n — 1 or 2, the result is easily verified. For n > 3, we construct for 
each value of n (mod 4) a suitable permutation. Figure [2] shows an example 
for each case of the proof. 

1. if n = 3 (mod 4), then any 2-permutation ir in S n is a valid candidate: 
indeed, W contains in this case exactly {n + l)/2 cycles of length 2, and 
Theorem 15.21 yields 

P*d(7r) > r— ^- - 1 = = — — . 



2. if n = (mod 4), we build a permutation a in S n by inserting a new 
first element as a fixed point in W, where ir is the permutation in S n -i 
constructed in the previous case, a contains n/2 cycles of length 2 
and one cycle of length 1 that corresponds to the fact that o~i — 1. 
Theorem 15.21 then yields 

_ . , n + 1 + ^ + 1 2n + 2 + n + 2-4 3n 

ptd(a) > 1 = = — . 

v V ; ~ 2 4 4 

3. if n = 1 (mod 4), we build a permutation £ in 5 n by inserting a fixed 
point anywhere in a, where a is the permutation in S n -\ built in the 
previous case. £ contains (n + 1 — 2)/2 cycles of length 2 and two cycles 
of length 1, and £i = 1. Theorem 15.21 then yields 

. . n + 1 + + 2 2n + 2 + n + l- 2 + 4- 8 3n-3 

ptrf(0 > 2 = = . 

F - 2 4 4 

4. if n = 2 (mod 4), we build a permutation r in by appending a 3- 
cycle to any permutation 7f such that 7r is a 2-permutation in S n ^%. f 
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contains (n + 1 — 3)/2 cycles of length 2 and one cycle of length 3, and 
Theorem 15.21 yields 




n = 3 (mod 4) n = (mod 4) n = 1 (mod 4) = 2 (mod 4) 



Figure 2: Cycle graphs of permutations it in S n with ptd(ir) > |_3n/4j, for 
all values of n mod 4. 

We can actually show that the lower bound on the prefix transposi- 
tion distance of 2-permutations is tight. In order to do that, we will need 
the following result. We use the following relation to order black arcs: 

(7Ti,7ri_i) -< (TTjjTTj-l) if j > L 

Lemma 5.5. jl/ For any tc in S n , let C\ be a cycle of length 2 in G(tt) with 
black arcs a, b; then there exists another cycle C 2 in G(n) containing two 
black arcs c and d such that a~<c~<b~<dorc-<a~<d~<b. 

This result can be interpreted in a more visual way in the case of a 2- 
permutation it by saying that in G(ir), every 2-cycle intersects with another 
2-cycle. We are now ready to prove the following result. 

Proposition 5.1. For any 2-permutation tt in S n , we have ptd(ir) = (3n — 
l)/4. 

Proof. The lower bound has already been observed in Theorem 15.41 To 
show that it is also an upper bound, we give an algorithm that sorts tt in 
exactly that number of steps. By Lemma I5.5[ every 2-cycle intersects with 
another 2-cycle, and as observed by Bafna and Pevzner 0], a sequence of 
two transpositions on any two crossing 2-cycles will transform them into four 
adjacencies: 
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7Tl 



7T,_l 7To 



becomes 




7T, 



7Tj_l 



7To 



which becomes 



7Tj 7T fc „! TTi TTj-l TTl 7Tj_i 7Tfc 

We transform the leftmost 2-cycle and any 2-cycle it crosses into four 
adjacencies using two prefix transpositions, which transforms n into a per- 
mutation a that contains r ^ — 2 cycles of length 2 and fixes the first element. 
Then, we carry out again this process until a is sorted, but we need three 
prefix transpositions at each step, since one move must be wasted to move 
the fixed points in cr's prefix out of the way, for instance as follows: 




The algorithm is guaranteed to terminate, since after applying each sequence 
of three transpositions of the form described above, we obtain either t, or a 
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permutation on which we can repeat the same process by Lemma [5.51 The 
proof follows from the fact that the number of prefix transpositions used by 
this algorithm is 

3 fn+1 \ 8 + 3n-9 3n-l 

2 + 



2 V 2 / 4 4 

□ 



6 Experimental results 



We generated all permutations in S n , for 1 < n < 12, along with their prefix 
transposition distance, and compared lower bounds ([3]), (ED and ([5]) to the 
actual distance. Table Q] shows the results. It can be observed that many 
more permutat ions are tight with r espect to our l ower bound (column 5) tha n 
with respect to Dias and Meidanisl ' (column 3) or Chitturi and Sudborough 's 
(column 4). 



n 


n\ 


tight w.r.t. (ED 


tight w.r.t. 05) 


tight w.r.t. (ED 


1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


3 


6 


4 


4 


6 


4 


24 


13 


15 


22 


5 


120 


41 


48 


106 


6 


720 


196 


255 


574 


7 


5 040 


862 


1 144 


3 782 


8 


40 320 


5 489 


7737 


27471 


9 


362 880 


31033 


44187 


229167 


10 


3 628 800 


247006 


369 979 


2103 510 


11 


39 916 800 


1706 816 


2 575 693 


21280 564 


12 


479 001600 


16 302 397 


25 791862 


236 651919 



Table 1: Comparison of all known lower bounds on the prefix transpositio n 
distance. Column 3 lists the number of cases where Dias and Meidanis's 



lower bound is tight 13l. p age 48] , column 4 lists the number of cases where 
Chitturi and Sudboroughr s lower bound is tight, and column 5 lists the num- 
ber of cases where our lower bound is tight. 

We also examined how large the gap between our lower bound and the 
actual prefix transposition distance can get. Table [2] counts permutations 
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whose prefix transposition distance equals our lower bound plus A. We note 
that, for n < 9, all permutations have a prefix transposition distance that is 
at most our lower bound plus 2 (plus 3 for n < 12). 



Tl 




A - 


A = 1 


A = 2 


A - 3 


l 


1 


1 











2 


2 


2 











3 


6 


6 











4 


24 


22 


2 








5 


120 


106 


14 








6 


720 


574 


143 


3 





7 


5 040 


3 782 


1234 


24 





8 


40 320 


27471 


12 310 


539 





9 


362 880 


229167 


128 576 


5137 





10 


3 628 800 


2103 510 


1427966 


97321 


3 


11 


39 916 800 


21280 564 


17532 948 


1 103 254 


34 


12 


479 001600 


236 651919 


221680 237 


20 667140 


2 304 



Table 2: Number of cases where our lower bound underestimates ptd(n) by 
A. 



7 Further observations on tt 

Now that we have an alternate representation of the cycle graph of a per- 
mutation as another permutation, we would like to examine whether or not 
other results can be obtained that could be helpful in getting insight on 
problems related to length-constrained factorisations of permutations. We 
investigate in this section a few relations between n and W, starting with 
relations between the cycle structures of both permutations when subjected 
to particular operations. 

7.1 Cycle structures 

A natural question is whether conjugacy classes are preserved by /(•), i.e. 
whether W and 7r cr are in the same conjugacy class for any choice of 7r and o 
in S n . The answer is negative in general, as the following counter-example 
shows: 7r = (1,3)(2) and r = (1,2) (3) are conjugate, but W = (0,1,2,3) o 
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(0, 1, 2, 3) = (0, 2)(1, 3) and r = (0, 1, 2, 3) o (0, 3, 1, 2) = (0)(3, 2, 1) are not. 
However, the relation we are interested in holds for two particular cases, 
whose significance we explain below. 



The following result is similar in spirit to iTannier and Sagotrs characteri- 
sation of "inverse breakpoint graphs" of signed permutations 24] . and shows 
that the cycle graphs of a permutation and of its inverse have exactly the 
same cycle structure. 

Lemma 7.1. For any tt in S n , we have tt~ 1 = (rf^ 1 )^ 1 \ 
Proof. Straightforward: 

~ = (0,l,2,...,n) o^Tr-STr-V--,^ 1 ) 



TT 



TT 1 O 7T O (0, 1, 2, . . . , n) O 7T 1 O (0, 72, U — 1, . . . , 1) O TT 
7T _1 O (0, 7Ti, 7T 2 , . . . , 7T n ) O (0, U, U ~ 1, . . . , 1) O 7T 



-1 1 

7T O 7T O 7T 



(Tf- 1 )^ 1 ). 

□ 



Tannier and Sagotr s idea of examining how the cycle graph of 7r _1 evolves 
when applying a signed reversal to a permutation tt - which reverses and flips 
the signs of the elements of an interval of it - was a key point in their suc- 
cessful attempt at designing an algorithm with an improved running time 
for sorting permutations by signed reversals. The above relation allows us 
to derive a simple description of the more general situation (albeit restricted 
to "traditional", unsigned permutations), i.e. how 7r _1 changes when an ar- 
bitrary rearrangement a is applied to tt. 



Corollary 7.1. For all tt, o in S n , we have (tt o a) 1 = (a 1 o tt 1 ) cr 
Proof. Lemma [3.11 yields 



(tt o a) 1 = a 1 o 7r 1 = (j 'oi 



a oct octoct o7t 1 oct (using Lemma 17.11) 



(If- 1 O TT" 1 )"" 



□ 
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A second particular case of conjugate permutations whose transformation 
by /(•) yields two conjugate permutations is presented below. Recall that x 
is the reverse permutation, i.e. \ = (n n — 1 • ■ ■ 1). 

Observation 7.1. For any 7r in S n , we have ix x = ((7f- 1 )x)(°> 1 ' 2 >--->' 1 ). 
Proof. We have by definition: 

W = (0,1,2,. ..,n)o(0,7r* tt*^,...,^) 

= (0, 1, 2, . . . , n) o n x o (0, n, n - 1, . . . , 1) o (vr*)" 1 

= (tt* o (0, n, n - 1, . . . , 1) o (tt*)- 1 o (0, 1, 2, ... , 

= ( X o tt o (0, 1, 2, . . . , n) o vr" 1 o x o (0, 1, 2, . . . , n)^ ' 1 ' 2 --") 

= ( X o vr o (0, 1, 2, . . . , n) o vr" 1 o (0, n, n - 1, . . . , 1) o x f^,...,n) 

= (xor^xp---) 

□ 

Conjugating vr by x corresponds to computing its reverse complement: in- 
deed, TTOX = (ll n 7T n _i • • • 7Ti), and Xo(vrox) = (n + l-7T n n+l - 7T n _! • • • n + 1 - 7Ti). 

By definition, 7r and 7r x have the same cycle structure, and by the above 
result, so do their images by /(■). The reverse complement operation is 
interesting because most (but not all, prefix distances being notable excep- 
tions j2l[) genome rearrangement distances are, in addition to being left- 
invariant, also "reverse complement-invariant": for all 7r and a in S n , we 
have d(ir, a) = d(Tc x , a x ). As a consequence, bounds obtained on the distance 
between 7r and i with respect to a certain set of operations can sometimes 
be improved by examining 7r _1 or tt x . 



Eriksson et al. introduced another important equivalence relation on 
permutations that does not preserve their cycle structure in the classical 
sense but that does preserve the cycle structure of their cycle graphs. This 
equivalence relation, whose equivalence classes are called toric permutations, 
proved useful in improving bounds on the transposition distance |ll|, M. We 



will see below that W provides a simple way of navigating through all cycle 
graphs of the permutations in the same equivalence class. The equivalence 
relation uses the following notion. 

Definition 7.1. The circular permutation obtained from a permutation n 
in S n is 7T° = 7i"i 7T2 • ■ ■ 7r n , with indices taken modulo n + 1 so that 

= < = 
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This circular permutation can be read starting from any position, and the 
original "linear" permutation is reconstructed by taking the element following 
as 7Ti and removing 0. For x in {0, 1, 2, . . . , n}, let x" 1, — (x+m) (mod n+1), 
and define the following operation on circular permutations: 

™ i „o T\ m -=— m -=— m 

111 + 71 = U 7Ti 7T2 ■ ■ ■ 7l n . 

Definition 7.2. For any tt in S n , the toric permutation tc° is the set of 
permutations in S n reconstructed from all circular permutations m + tc° with 
< m < n. 

Definition 7.3. Two permutations ir, a in 5*^ are torically equivalent if 
a G 7i° (or 7r G cr°), which we also write as tt =° o. 

Let us illustrate those notions using our running example n = (416257 3); 
we have tt° = 4 1 6 2 5 7 3, and 



+ 7T° 


= 04162573 


1 +7T° 


= 15273604 


2 + 7T° 


= 26304715 


3 + 7T° 


= 37415026 


4 + 7T° 


= 40526137 


5 + 7T° 


= 51637240 


6 + 7T° 


= 62740351 


7 + 7T° 


= 73051462 



which yields < = {(4 1 6 2 5 7 3), (4 1 5 2 7 3 6), (4 7 1 5 2 6 3), (2 6 3 7 4 1 5), 
(5 2 6 1 3 7 4), (5 1 6 3 7 2 4), (3 5 1 6 2 7 4), (5 1 4 6 2 7 3)}. Hultman 



171 ] proved the following interesting result. 



Lemma 7.2. fil l For all it in S n and < m < n: every cycle in G(tt) 
is mapped onto a cycle in G(a), where a is the permutation obtained from 
m + 7r°. 

In other words, if it =° a, then W and a are conjugate. We show below 
how one can iterate over the cycle graphs of all elements in tt°. 

Lemma 7.3. For all it, a in S n : if a° = m + tt° , then a = W (°> 1 > 2 >-> n ) m . 
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Proof. By Equation [IJ we have: 

a = (0, 1,2, . . . ,n) o (a ,cr n ,a n __i, . . . ,a x ) 

= (0, 1, 2, . . . , n) o (m + vr , m + 7r n , m + 7r n -i, • • • , m + 7Ti), 

since by hypothesis cr° = m + 7r°. 

On the other hand, the mapping (ttq, ir n , n n -i, . . . , tt\) !->■ (m + 7r , m + 
7r n , m + 7r n _i, . . . , m + 7Ti) consists in replacing each element of the cycle with 
its value plus m (mod n+1), which is by definition equivalent to conjugating 

(7T , 7T n , 7T n _i, . . . , 7Ti) by (0, 1,2,..., Tl) m . □ 

Corollary 7.2. For a// vr in S n , we have {a \ a e <} = {tH ' 1 ' 2 --")" 1 | < 
m < n}. 

Other relations between the cycle structure of 7r and that of W can easily 
be derived from previous work. The following relation allows us to bound 
the number of odd cycles of W. 

Theorem 7.1. [20] For all it in S n , we have td(ir) < n — c 0( id(X(ix)) . 

The following result is an immediate corollary of Theorems 14.21 and 17.11 

Corollary 7.3. For all tc in S n , we have 2c oc y(r(7r)) < n — 1 + c ^(r(7f)). 

Similarly, the following result is an immediate corollary of Theorem 14.11 
and of the characterisation of exchanges as restricted block-interchanges. 

Corollary 7.4. For all n in S n , we have 2c(T(n)) < n - 1 + c(T(n)). 
7.2 Descents of it and cycles of W 

Aside from relations between cycle structures, we can also establish relations 
between pairs of elements of 7r and cycles of W. An example of such a relation 
is the fact that the number of adjacencies in (0 7Ti 1x2 • • • ir n n + 1) equals 
ci(r(7f)). We will prove that a less obvious relation connects the descents of 
7r (defined below) and the cycles of W. 

Definition 7.4. A descent in a permutation it is a pair (7Tj_i,7rj) such that 

TTj < 7Ti_i. 

For instance, the permutation (4|1642574 3) has three descents, 
indicated by vertical arrows. 
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Definition 7.5. A cycle C in G(tt) contains a descent (7Tj_i,7rj) if (71$, 7Ti_i) 
is a black arc of C . 

We now derive bounds on the number of descents contained by cycles in 
G(tt). 

Lemma 7.4. For all n in S n , every cycle of length £ > 2 in G(tt) contains 
at most £ — 1 descents and at least one descent of it. 

Proof. For clarity, let us write the vertices of C in the order in which C 
visits them, starting with the element whose position in ir is maximal: we 
get C = (7^,7^, 7Tj 2 , it j2 , . . . , n ik ,7i jk ), where it (resp. j k ) is the largest (resp. 
smallest) position of an element of tt appearing in C. We identify here tt 
and 7r n+ i = 7T0 (mod ri + 1). Recall that (n ix , iTj x ) for 1 < x < k is a black 
arc of C, and that by Definition I2.10[ the following relation holds: 

7Ti x = 7Tj :r _ 1 + 1 for all 1 < x < k, and = 7Tj k + 1. (6) 

1. for the upper bound: assume on the contrary that C contains £ de- 
scents; then every black edge of C corresponds to a descent, and we 
have: 

TT jx > ir ix for 1 < x < k. (7) 

By alternating between the conditions specified by Equations [7] and El 
we obtain: 

*3k > n ik = ^'fc-i + 1 > ^ifc-i + 1 = Ti*-a + 2 > • • • = TTfo + k — 1 

> -K ix + k — 1 = 7rj fc + fc, 
which is clearly a contradiction. 

2. for the lower bound: assume on the contrary that C contains no de- 
scent; we have: 

TTj x < 7i ix for 1 < x < k. (8) 

By alternating between the conditions specified by Equations and El 
we obtain: 

^1 - 1 = Kj h < *ik = ^ifc-i + 1 < Kifc-x + 1 = 7T jfe _ 2 + 2 < • ■ ■ = 7T h + k ~ 1 

< n h +k-l. 
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For the above relations to hold, elements from the set A = {vr ife , 7Tj fc _ 1 + 
1, . . . , 7Tj 2 + k — 2} can only be assigned values from the set B = {n^ + 
1, 71"^ +2, . . . , 71"^ + k — 2}. However, we have k — 1 = |A| > |£?| = k — 2, 
which clearly makes it impossible to obtain a permutation. 

Finally, note that it and 7r can be regarded as equivalent as far as descents 
are concerned, since (ttq,tti) and (7r n ,7r n+ i) cannot be descents. □ 

The following result is a direct corollary of the above. 

Proposition 7.1. For any 2-permutation tc in S n , we have des(n) — (n + 
l)/2. 

Proof. By definition, W contains exactly (n + l)/2 cycles of length 2, and by 
Lemma [7. 4[ each of these cycles contains exactly one descent of tt. □ 

8 Conclusions 

We presented a new framework for reformulating any edit distance problem 
on permutations as a minimum-length factorisation problem on a related 
even permutation, under the implicit assumption that the edit operations 
are revertible. This approach is based on a new representation of a structure 
known as the cycle graph, which pervades the field of genome rearrangements 
in several different forms; it previously allowed us to enumerate permutations 
whose cycle graph decomposes into a given number of alternating cycles (lo| . 
and allowed us in this work to recover two previously known results in a 
simple and unified way. Moreover, we used our approach to derive a new 
lower bound on the prefix transposition distance that, as we showed both 
theoretically and experimentally, is a significant improvement over previous 
results. From that result, we deduced an improved lower bound on the prefix 
transposition diameter of the symmetric group, whose exact value is still 
unknown. Finally, we investigated other relations between permutations and 
their cycle graphs that we hope will prove useful in obtaining new results. 

Several interesting questions and leads for future work arise. First, our 
method provides an automated way of obtaining lower bounds on distances 
between permutations; is there an analogous way of obtaining upper bounds 
instead? Second, we initiated the study of relations between statistics on a 
permutation and statistics on the permutation that corresponds to its cy- 
cle graph. Can other relations be deduced and used to prove other results, 
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including tighter bounds on the distances of interest? Third, permutations 
are but one structure for which the cycle graph has been defined. Other 
structures, such as signed permutations, give rise to a more general structure 
known as the breakpoint graph. Are there analogs, or generalisations of /(•) 
that can yield similar results on signed permutations as well? Finally, another 
question is whether Cayley graphs obtained from genome rearrangement op- 
erations can yield good interconnection networks. For instance, (signed) re- 
versals generalise the operations that generate the (burnt) pancake network, 
exchanges generalise the operations that generate the star network, and prefix 
transpositions generalise the operations that generate the bi-rotator graphs 
(see Lakshmivarahan et al. 23[ for definitions). It seems likely that collab- 
orations between researchers in both fields could be fruitful in investigating 
this topic. 
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